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ABSTRACT 


Two design methods for multi-parameter dynamic systems 
are proposed. They are intended to eliminate the limita- 
tions and disadvantages of the existing asian methods. The 
powerful mathematical tools of optimal control theory are 
applied to the practical design problems of classical con- 
ies 

The first method is intended for linear systems only; 
the design problem is solved in the s-domain, by finding 
"tHe best root locations" of the system's characteristic e- 
quation. In the second method, the -design problem is solved 
by finding "the best response" of the system in the time 
domain. The second method is applicable to a wide range of 
dynamic systems; it can be used to synthesize linear, non- 
linear and sampled-data systems, and systems with time de- 
lay. This method is also extended to a numerical stale, 
analysis procedure. 

Fourteen examples are presented to illustrate the 


applications of the methods. 
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I. INTRODUCTION 


Classical control theory and modern control theory 
have been extended in diverging routes since the inception 
of modern control theory. Modern techniques have been suc- 
cessful for larger scale problems and for problems which 
include some new performance criteria, such as expenditure 
Of Control energy, etc., but a general purpose or a special 
purpose computer in the control loop is usually necessary 
to implement the solution. 

Most design techniques of classical control theory 
were developed before the ot teieive application Of "areca. 
computers to system design problems. Using a digital com- 
puter in the design phase of a control system is usually 
feasible, especially if the design can be implemented by 
using passive circuit elements rather than complicated and 
expensive memory units. 

During the last decade mathematical methods and numer- 
ical techniques have been developed and applied to the prob- 
lems of "Optimal Control Theory". The same methods have 
found application in other areas of science, but they have 
not been applied efficiently to the design procedures of 
the dynamic systems which have been defined in the sense of 

classical control theory. The classical techniques, al- 
though having many drawbacks and limitations (see Section 
II.A), were well developed by the time effective optimiza- 


tion methods arrived on the scene. The classical methods 


have been successful for a restricted class of design prob- 
lems; for this reason probably, modern control theory grew 
independently, rather than improving and replacing the 
existing methods. 

In this thesis it is intended to narrow the gap be- 
tween the classical and modern control theories. Powerful 
mathematical tools of optimization theory are applied to 
the practical design problems of classical control. Two 
general design methods for multi-parameter dynamic systems 
are proposed: 

1. Optimization for the Best Root Location in the 
s-Domain, and 

2. Optimization for the Best Response in the Time 
Domain. 

Chapter II of the thesis is devoted to the first meth- 
od which is intended for linear systems only -- the design 
problem is solved in the s-domain by finding "the best root 
‘seer ere” of the system's characteristic equation. Four 
linear system design examples are included. The rest of 
the thesis is devoted to the second method and its various 
applications. In Chapter III, a general description and 
the philosophy of the method "Optimization for the Best Re- 
sponse in the Time-Domain" is given, and it is applied to 
the same linear system design problems discussed in Chapter 
II, to provide a basis for comparison of the results ob- 


tained by both methods. 


The time-domain method is applicable to a large varie- 
ty of system design problems. To demonstrate this point, 
Chapter IV is devoted to various applications of the method. 
Three examples for nonlinear and sampled-data systems and 
a system with transport lag are presented in Chapter IV. 

[nsGhnapter Vv, a numerical method to investigate tiie 
stability limit of the dynamic systems, is introduced as an 
extension of the time-domain method. 

In the year 1964, the gap between control theory and 
the practical control problems became substantial and ex- 
tensive research effort has been directed to bridge this 
gap between the theory and paractice (see, for example, Ref. 
2). The methods presented in this thesis can be thought of 
as another approach to suboptimal system design problems. 

In this respect, these methods can be considered as applica- 
tions and extensions of Optimal Control Theory; on the other 
hand, these methods provide better solutions to the practi- 
cal dynamic system synthesis problems of Classical Control 
and they can be considered as the extensions of Classical 
Control. Since the extensions from Optimal Control Theory 
and Classical Control are toward each other, these methods 
hopefully provide one of the necessary bridges to reduce 


the gap. 


II. GLINEAR SYSTEM DESIGN BY OPTIMIZATION FOR 


THE BEST ROOT LOCATIONS IN THE S-DOMAIN 


A. GENERAL 

Present design methods of classical control theory are 
limited and have some serious disadvantages. Root-locus 
and frequency-domain techniques, for example, can handle 
systems which have only one parameter. Algebraic methods 
(in the coefficient and parameter planes) are two-variable 
procedures. The major design tools of these methods are a 
set(s) of graphs, and obtaining these graphs for higher- 
order systems is quite a laborious task. The trial-and- 
error nature of the classical design procedures is, however, 
the most significant limitation of these methods. 

Since the dynamic behavior of second-order systems has 
been thoroughly studied and is well known, most of the pre- 
sent methods for high-order system design depend on the 
dominant-roots concepts. Finding a pair of complex dominant 
roots for a high-order system is equivalent to approximat- 
ing the system by a second-order system. 

A pair of complex roots of high-order system can be 
made dominant by satisfying at least one of the two condi- 
tions given below: 

1. Transients in the time response due to undesired 
roots (all the roots except the selected "dominant" pair) 
last for a much shorter time than the transients of Cheweas 


lected pair. This is achieved if the time constants of the 
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undesired roots are much smaller than those of the selected 
pair. In the s-domain, the magnitudes of the real parts of 
the undesired roots must be much larger than the magnitu- 
des of the real parts of the selected pair to eeion this 
condition. The remote undesired roots may cause undesir- 
able effects in the time response for a short period of 
time (near zero time), but thereafter the time response is 
governed by the selected roots. 

2. The coefficients in the time response associated 
with the desired roots are much larger than those associ- 
ated with the undesired roots. In other words, the contri- 
bution of the selected pair to the transient response is 
much larger than the contribution of the other roots. In 
the s-domain, this condition is satisfied if the residues 
of the desired roots are much larger than the residues of 
the other roots. 

If both of these conditions are satisfied, the time 
response of the selected pair dominates the overall time 
response, and the selected pair of roots is called 
"dominant". If only one of the conditions is satisfied an 
acceptable transient response may still be obtained. The 
first condition alone, for example, has been used in the 
s-domain design procedures to avoid the lengthy and time- 
consuming residue calculations at each step of the trial- 
and-error procedures. 

When designing a high-order feedback system by root- 


locus methods, or by algebraic methods, the designer 


iG 


usually selects a pair of complex root locations in the s- 
domain, and then tries to locate two of the roots of the 
system's characteristic equation at these selected loca- 
tions and the other (undesired) roots at remote points of 
the left-hand side of the s=plane. If this can be achieved, 
the transient response of the system closely follows the 
selected (desired) second-order system's response. There 
are usually other specifications to be met, such as the 
steady-state error, maximum allowable gain, and limitations 
on the system components, which are treated as additional 
constraints. 

The number of free independent parameters (No) cannot 
exceed the order of the characteristic equation (n), and if 
the number of free independent parameters is equal to the 
order of the system (NV. =n), the roots of the characteris- 
tic equation can be located at any point in the s-plane. 
The design is completed by simple algebraic manipulations 
(see,for example, Ref. 12). If the number of free indepen- 
dent parameters is less than the order of the characteris- 
tic equation of the system (NO <n), however, the designer 
loses his control over some of the roots. Since one free 
parameter exactly locates one root, when all parameters are 
used to locate an equal number of roots, the excess (n-N) 
root(s) may be anywhere in the s-plane. Some of these 
roots may move very close to the origin and become dominant, 
or some may even go into the right half of the s-plane and 


make the system unstable. 
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Two parameters can be used to locate a pair of complex 
roots exactly at the specified locations, but these roots 
can be constrained to move on a line by using just one pa- 
rameter. For example, their real parts (0) or imaginary 
parts (W) can be fixed, or they can be forced to move on a 
Constant damping ratio (€)-line or on a constant natural 
frequency (w)-circle. 

When designing with a fewer number of parameters than 
the degree of the characteristic equation, the procedure is 
a trial-and-error method and success may depend primarily 
on the past experience of the designer. There are several 
methods for locating two of the roots of the characteristic 
equation at the desired locations, but if the number of 
free parameters is less than the order of the system, none 
of these methods ensure 

a. the stability of the system, and 


loys the s@ominancesoft the selected» roots. 


Be PHILOSOPHY OF OPTIMIZATION 

Since placing a pair of complex roots at the desired 
locations has no significance without ensuring their domi- 
nance and the system's stability, the method proposed in 
the first part of the thesis provides these two conditions 
by confining the undesired roots to appropriate region(s) 
of the s-plane, and then finds the optimum parameter values 


to locate the dominant roots as close to the specified 


iss 


locations as possible’. Since the undesired roots are not 
placed at specific locations (they are only confined to 
stay in a remote region of the s-plane) none of the free 
parameters are "used up"; only their freedom is limited. 
The problem is solved in the "Parameter Space" which 
is a N,-dimensional Fuclidean space. A performance index 
(J), which is usually a measure of closeness of the dominant 
roots to the desired root locations, is defined. By Gee 
mizing the performance index with respect to the selected 
free parameters (z) (see Section II.D) one obtains the 
"best possible solution" to the design problem as it is 
specified. The design specifications yield some linear and 
nonlinear constraining equations and/or inequalities; it 
is assumed that these constraining equations define a bound- 
ed region in the N.-dimensional parameter space. Some of 
the constraining equations are generally nonlinear and it 
may sometimes be impossible to find a feasible point (a 
point which simultaneously satisfies all of the constrain- 
ing equations and inequalities) to minimize the performance 
index. In practical terms, this means that the design (or 


compensation) is not always attainable by using a certain 


t For the simplicity of explanation it is assumed here 
that the design is carried out for a selected pair of 
complex roots; but this is not mandatory. The method is 
generalized (in Section II.B) and an example is given 
(Example 3) for the design of systems without a pair of 
selected dominant roots. 


type of compensator. If this occurs, the designer should 
change the system configuration (usually by changing the 
type of compensator, or by inserting additional compensa- 


Eilon) . 


Ce CONSTRAINTS 

Performance specifications and available system com- 
ponents dictate the constraints, which can be divided into 
three general groups: 

1. Constraints due to physical limitations, such as 
limitations on the available system components. 

2. Constraints due to static (or steady-state) per- 
formance specifications. 

3. Constraints due to dynamic performance specifica- 
tions. 

The constraints in the first two groups are linear or 
can be approximated by linear inequalities. The constraints 
in the last group are obtained by factoring the system's 
characteristic equation (as explained in the next section) 
and they are expressed as "n" linear and/or nonlinear si- 


multaneous algebraic equations. 


Dis FACTORING THE CHARACTERISTIC EQUATION 

The denominator of the system's transmission function 
(from input to output) is called "the characteristic poly- 
nomial", and when it is equated to zero it is called "the 


characteristic equation". The characteristic polynomial of 


£5 


the system is a n)_degree polynomial in the complex vari- 
able "s", and a function of N. free system parameters (Q@) ; 


that is: 
ee 
_ al 
Pon's) 7 By 2 , (Ii ow 
=0 


&a& s 8 e 
where a = 1 and other coefficients are either known con- 
stants or functions of one or more system parameters. 


Therefore, the characteristic polynomial can be represented 


as 
Pan(S) = Popes) - (II.2) 


To represent the desired and undesired roots explicit- 
ly, it is necessary to factor the characteristic polynomial 
into two parts. In this factored form the characteristic 


polynomial may be written as 


Py (& 8) = | P(o,w,s) | ° |P,4(B.s) | (-TDaiews3 ) 


where f represents the coefficients in Pid (a polynomial of 
degree "n-2" ins, representing the undesirable roots) and 
Py 1S a quadratic in s which represents the desired roots. 


Two forms of P, will be considered: 


“~ 


P,(0,W,s) = s° + 208 +(0* + uw"), (IT:4) 
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where 0 and W are the real and imaginary parts of the 


desired complex roots, and 


es 2 
P4(o,%,Ss) = a 20 Ww Ss + Ww F (Tres) 


where € is the damping ratio and Ww 1s the natural frequen- 
cy of the second-order system. 

The characteristic equation can be factored by using 
several methods, such as, by dividing the characteristic 
equation by Py Ox Pid and equating the remainder (in each 
case) to zero. In another method the general forms of Ps 
and a are multiplied together and the coefficients of the 
product polynomial are equated to the corresponding coef- 


ficients of the characteristic polynomial; that is, the 


identity shown below is solved 
Pa - PSa = Pan ; (IT.6) 


These processes yield n algebraic equations, all of 
which must be satisfied to obtain the desired factored 
form of the characteristic equation. 

The left side of Eq. (II.3) contains N. free system 
parameters (@), and the right side is a function of n root 
parameters (0, 8); that is, there are a total of N. + n 
parameters available, but only No of them are independent -- 
the others are related to the independent ones through the 
n algebraic constraining equations which are obtained during 


the factoring process. Any, N. of these parameters can be 


IL 


chosen as free variables (Z); the remaining n variables 
are determined by the constraining equations due to dynamic 


performance specifications (see Paragraph II.C(3) above). 


BE. THE PERFORMANCE INDEX (J) 

When the characteristic equation is factored, that is, 
after the free variables (Z) are chosen and the constrain- 
ing equations are found, a performance index (J) is defined 
so that when it is minimized with respect to the chosen 
variables (Zz) in the bounded region defined by the con- 
straining equations and inequalities, it yields "the best 
root locations". The best root locations may mean that the 
desired roots are as close to the specified locations as 
possible (if the dominant roots concept is being used in 
the design), or to satisfy other criteria -- as explained 
below -- and that the undesired roots are in remote regions 
of the left-hand side of the s-plane (see Figs. II.1 and 
Tis237 

The form of the performance index depends on the de- 
sign specifications. If the design is to be carried out by 
locating a pair of dominant roots as close to the desired 


locations as possible, then the performance index may take 


the form 
a (O 5 ~ o)* + (ww - w)* - ao (Tima 
T= (0, -0)% + (w - w)? (II.8) 
d 


where the subscript "d" indicates the desired coordinates. 
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Fig. II.1. The desired and undesired root location in 
the s-plane. 


NOTES: Point (A) is an actual root location at a 
minimization step (J = ry ) 
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Fig. II1.2 The desired and undesired root locations in 
the s~plane., 


Sa-d, @z-.¢0 


(Point (A) is an actual root location ata 
minimization step). 


Ike 


Sometimes design specifications allow the designer to 
use a simpler performance index. For example, only the 
first terms of Eq. (I1.7) or (II.8) may be used. For some 
purposes extra terms may be added to the performance index. 
For example, it may be desired to confine the destred@eege. 
to two small regions around the desired root locations and 
the real part of the closest undesired root (to the origin) 


is maximized by minimizing the performance index. 
J =i ’ (q ? 0) (THs9) 


where q is the magnitude of the real part of the closest 
undesired root. The negative sign in (11.9) converto game 
maximization problem into a minimization problem. 

‘For different kinds of dynamic performance specifica- 
tions, different performance indices can be defined and 
minimized with respect to the chosen free parameters. If 
the dynamic performance specifications are very tight (that 
is certain rise time, overshoot, settling time, frequency 
of the transients, etc., are desired) it is best to use the 
dominant-roots concept (Eq. (1.7) or (I1I.8)) for the per- 


formance index. 


I PROGRAMMING 

Although it is possible to solve simple problems by 
hand calculations, in general, a computer program capable 
of minimizing a nonlinear function with linear and nonlinear 


constraints is necessary for higher-order systems. For the 
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examples considered in this part of the thesis SUBROUTINE 
BOXPLX (see Refs. 1 and 4) has been used in the minimiza- 
tion with good results. A computer program which is written 
in FORTRAN IV, and used for the minimization part of Ex- 
ample - 4 is given at the end of the thesis (see COMPUTER 


PROGRAM 1). 


Ge EXAMPLE lL 
A linear third-order feedback control system is shown 
in block diagram form in Fig. II.3(a). A step-by-step pro- 
cedure to compensate the system for the desired transient 
response by using the method "Optimization for the Best 
Root Locations in the s-Domain" is given below. 
1. Interpret the Desired Dynamic Performance Specifi- 
cations 
The desired dynamic performance specifications are 
usually given as maximum overshoot, rise time, settling time, 
frequency of the transients, etc., and it is usually best 
to find a pair of complex root locations in the s-domain to 
represent these specifications. For this problem it is as- 


sumed that a pair of caqmplex roots located at 0, = 0.5 and 


d 


Wa = + 0.5 would yield the desired dynamic response. 


2. Check to Determine if the Desired Roots are on the 
Root Locus of the Uncompensated System 
The designer gains considerable insight into the 
design Bromlen by just sketching the root locus; this also 


helps in deciding on the type of compensation that should 


ZA 


be used, if different choices are available. The root 
locus for the uncompensated system is shown in Fig. I1I.3(b), 
and the desired locations are not on the root locus. 
3. Decide on the Type of Compensation 
A cascade compensator will be used. The system 
with the cascade compensator is shown in Fig. I1.4(a). 
4. Find the Transmission Function and the Character- 
istic Equation 
The compensated system becomes fourth order (n=4), 
and there are three free system parameters (NO = 3) -- the 
pole location (p), the compensation ratio (k) of the com- 
pensator, and the system gain (K). The vector representing 


the free system parameters is 


a= (k pK). (II. 10) 


Cav] 


The transmission function of the compensated system is 


T(s) _ CCS) 9 9 CRA sii @ur 


e(sity aero) + (K/k) (s+kp) 


and the characteristic equation is 


S° + (286)s° + (142p)6* 4 (okR epee cen on 


5 Factor the Characteristic Equgitgon 
Using one of the methods described in Section II.D 
factor the characteristic equation and find: 


d 


a. The fae a = 2” -degree polynomial (Pg) to 


represent the undesired roots. 


Z2 


jw 









(a) System Desired root lLocatio 
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| 
| 
| 
| 
| 
(SEE DED AE) ma) 
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(b) Uncompensated system's root locus 
Fig. II.3. EXAMPLE 1 
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is 
a) Compensated syste 
(a) p ystem ———— aos 


eee 
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(b) Root locus for the compensated system 


Fig. II.4. EXAMPLE 1 
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speoeiaie 
Pig (fs) Se eee eS (Il. len 


where Pp = (a Sie : 


b. A quadratic in s (P.) to represent the desired 


roots (the roots which will be made dominant by optimiza- 


Elon). 
2 2 2 
P ,(o,W,s) = NG ee Cue Cement cy) ae) (TI.14) 
c. The type -3 constraining equations by using 
Ie un (ANA (5)) 
Pp (@s) = |P,(o,,s) | - [Pp 4(Bs) | . (1.6) 


For notational convenience let c & (07 4s co 


= + (2+p)s> + (1+2p) s° + ({(ptK/k)s + Kp = 
( Tigl>) 


se + 2(atc)s> + (ena eet 2(actbo)s + be. 


Fquate the coefficients of the polynomials of the. above 


identity and obtain the n(=4) constraining equations 


p + K/K =m lac + bo) ( Tage ) 
Kp = be | (ig? ) 
1+ 2p = 4a0 +b+tc (Tes) 


2 ape = 2a + a) ( Tier) 


Therer arent otal weds N. +n (= 7) variables (@,K, imams 


a,W), and any N. (= 3) of them can be chosen as the free 
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variables to constitute the z-vector of the minimization; 
the remaining n (= 4) will be related to the free variables 
by the constraining equations. After selecting the free 
variables, the constraining equations should be rearranged 
So. Ehat: 

(1) the left side of each equation is one of 
the dependent variables, and 

(2) the right side of each equation must con- 
tain only the free variables and the dependent variables 
defined in the preceding equations. 

To avoid substitution, 9, a, and c are chosen as the 


free variables, hence 
- 
zZ= (a> Oo (e)r—s, (TiS20) 


and the constraining equations are solved for the remaining 


four variables 


p= 2 (a +o - 1) tsi 410), 
b=1+ 2p-c - 4ao (II.22) 
K = be/p Bint, Oe) 
k = K/(2 (ac + bo) - p) , (II.24) 


Ge sede ENemn ced nme MclmCONs tia dmtss 
a. The compensator will contain only passive ele- 


ments; hence, 


pee eo : (II.25) 
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b. A physical constraint is assumed on the sys- 


tem gain 


0.4 =K © 0.5 ( Tia2e) 


c. The compensation ratio of a cascade compensa- 
tor should not be greater than ten. 


O.1 = k s 10 (ILD yay 


d. To ensure complex dominant roots, 


Woe (Tae 
e. To confine the undesired roots to the lefeten 
s = aches = ~1.5 line, 
2 L. 
aU fer Sy oe des (IT.29) 


f. All free and dependent parameters must be pos- 


itive: hence 


ole an. O) (IT .30) 


7. Define a Performance Index 
For this problem the performance index defined by 


Fg. (II.7) is suitable, 


J = (¢. = a)? + (Ww. _ ws) - (II .7) 


d 
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8. Minimize the Performance Index 
Find the best root locations and the corresponding 
optimum system parameters (@*) by minimizing the performance 
index with respect to the chosen free parameters Zz = (a ac) 
in the bounded parameter space defined by the constraining 
inequalities (II.25) - (I1I.30) and the constraining Eqs. 
ait) — CL wo) . 
9. Relax the Necessary Boundaries if Possible 
Tf at the end of a minimization process the select- 
ed roots cannot be located in a close neighborhood of the 
desired locations (determined by a preselected norm -- see 
Eq. (I1I.31) below), and the result is a constrained minimum 
(that is, at least one of the variables is on the limiting 
boundary) then the limiting boundary(s) should be relaxed 
(if possible) and the minimization should be repeated. In 


this example the value was the only boundary which 


Tit 
limited the minimization process. The boundary on the un- 
desired roots was relaxed toward the origin by 0.1 steps 


Werte a 1. 


/ 2 er 107" Gir oT) 


was obtained. After each relaxation, SUBROUTINE BOXPLX was 
called and the performance index was minimized again to ob- 
tain optimum parameter values for the new boundaries. The 
results obtained at the end of each minimization are tabu- 


lated in Table II.1. Relaxation was stopped when the 


VF 


undesired root (#1) was no longer on the o boundary (see 


Pe 
last line of Table II.1); at the same step the stopping 
criterion on the performance index (Eq. (I1I.31)) was also 
satisfied. 

The repeated minimization process and boundary relax- 
ation were achieved automatically in the main program (see 
the computational flow chart -- Fig. II.5). 

10. Interpret the Results 


At the first minimization step the undesired roots 


were kept to the left of s = -0 = -1.5 line. The desired 


TE 
roots were probably dominant (since the transient due to 


th of the time 


EHoeundeSiredmecal woot sett lose déwn inwWyc 
of the transients caused by complex roots), but the domi- 
nant roots were not close to the desired locations and the 
transient response would not be satisfactory. At the last 
step of the relaxation the desired roots were exactly locat- 
ed but their dominance was questionable since an undesired 
real root was located close to the desired roots. 

To make the final decision on the parameter values, 
residues can be calculated; this would certainly remove 
the doubts about the dominance of the desired roots, but 
Since a digital computer is assumed to be available and 
digital simulation of the system avoids the lengthy residue 
calculations, it is best to obtain the time response of 
the system and compare it with the desired second-order 


response. 
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OUTPUT 
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Fig.II.7. EXAMPLE 1. System's time response (A) with 
optimum parameter values and 
desired response (B) to a unit 
step input. 


SZ 


Ec ttso 1S the response (to a Unit step input) for 
the first line of Table II.1; the desired second-order 
response is also shown. Fig. I1.7 shows the time response 
of the system (to a unit step input) for the parameter 
values tabulated as the last line of Table II.1, where the 
complex roots were exactly located at the desired locations 
co = 0, = o.5 and & = Ws = + 0.5) but a real root was close 
to the desired roots (at s = -0.94). Since the compensated 
system's response is very close to the desired response, 


the optimum parameter values tabulated as the last line of 


Table II.1 yield a very satisfactory design. 


H. EXAMPLE*2~ 
The linear fourth-order system shown in Fig. I1I8(a) is 
given. Steady-state performance specifications require the 


system gain to be equal to 


K 





(a) The uncompensated system (b) The compensated system 


Fig. II.8. EXAMPLE 2 


(or greater thah) 1000 (K 2 1000), but the root-locus 


study shows that the system is unstable if K > 193.23. 
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1. The dynamic performance speci fications= can be met 
by locating a pair of complex (dominant) roots at Coa in 


and W 5 = + 2, 

2. The desired root locations are not on the uncom- 
pensated system's root locus, and the system gain is much 
less than the desired value (for the part of the root locus 
in the left hand side of the s-plane); therefore, some kind 
of compensation is Wee asce 

3. A combination of tachometer and acceleration feed- 
back will be used to compensate the system. The configura- 
tion of the compensated system is shown in Fig. II8(b). 

4, The transmission function of the compensated system 
is 


SS = : P) (Pieoe) 


aoe 16s” +(65+KK_) s“+(50+KK, ) s+K 


WIA 
alu 


IGS) = 


and the characteristic equation is 
4 3 2 = 
s + 16s~+ (65+KK_)s + (50+KK,)s tk =U (ELS) 
5. The desired-roois quadratic is 


P (os, oe Ss) = a7 + 20s tec (FT 234) 


where 


Paragraph numbers indicate the corresponding steps in 
EXAMPLE 1. 
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The undesired-roots polynomial is 
2 
Pig (Be s) = s° + 2as +b (ie 3 5)) 


where 


B= (a b) 7 : 


Eq. (II.6) is used to factor the system's character- 


istic equation, 


A 3 


s + 16s~+ (65+KK_) s“+ (50+KK,)s + K = 


(iisee) 


see Dilan (bt cudaoee 2(bct+ac)s + be. 


By equating the coefficients of the corresponding 


Ecmms Of Ldentity (11.36), 


a+o:= 8 
b+ca+t 4a0 = 65 + KK, 

Lit. 39) 
2(bo + ac) = 50 + KK 


bc = K 


The z-vector, representing the independent variables 


in the minimization process, is chosen to be 


z= (gcc eS) 2 : ; (II.38) 


The constraining equations are 
a = 8 -0 
b = WW 

(JEaE 5 SiS) 
K_= (b + c + 4ao —- 65) /K 
K 


= (2.(bo + ac) == s0y7K 
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6. All the variables must be positive (ouszero), and 


the other constraints are 


K 2 1LO0O 
Os Kee. Oo 
a 
Os Ky = 10 (II.40) 


g* - om <0 


1. 
a- cee 2 Ory iG (II.40a) 
7. The performance index is 
2 
J” =" (1220) + 02 .-w) (Toned 1 ) 


8. SUBROUTINE BOXPLX was again used for the minimiza- 
ee Kel aee | 

9. None of the variables was on the limiting bound— 
aries; therefore, no relaxation was necessary. The result 
of the minimization is tabulated in Table II.2. 

10. The compensated system's response (to a unit step 


input) and the desired response are shown in Fig. II.9. 





LIMITING! UNDESIRED ROOTS | DESIRED ROOTS SYSTEM 
“ VALUE PARAMETERS 
FOR THE 
UNDE- 
SITRED 
(oO. } | Real Imag. | 
= part part ae Kr 
7 
ae 999 ; -199.585 ; -1.001! 42.002 taeoe 1676 | ~4197 
TABLE II.2 
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Fig.II.9. EXAMPLE 2. Time response (A) of optimally com- 
pensated system and the desired re- 
sponse (B) (to a unit step input). 


oO” 


ds EXAMPLE 3 

The linear system given in Example 2 (and shorn Io. 
Fig, I1.8(a)) will be compensated for the fastest response 
(the shortest rise time and the shortest settling time). 
The steady-state requirements and the type of compensation 
to be used are exactly the same as in Example 2. 

Since the only dynamic system specification is "obtain 
the fastest response", the actual form of the transient re- 
sponse and the location of the dominating roots (whether 
they exist or not) are of no interest to the designer. This 
example shows that the dominant-root idea, although vate 
useful for many design problems, is not mandatory for the 
system désign procedures. 

As mentioned earlier (in Section II.A) the s-domain 
design methods try to relocate the roots of the system's 
characteristic equation to obtain a suitable transient re- 
sponse. Most of these methods ignore the magnitude effects 
(to avoid the residue calculations), and they try to clas- 
sify the transient response as "the fast-varying part", and 
"the slowly-varying part". This procedure depends on the 
assumption that, "magnitude-wise all of the roots have the 
same importance in the transient response”. 

With the same assumption one may try to obtain "the 
fastest response" by making the real parts of all the 


roots as large as possible, because the rise time and the 


3 


settling time of a transient response caused by a root is 
inversely proportional to the real part of that root. 

The solution of the problem contains exactly the same 
steps as the previous example (see Eqs. (II.32) - (I1.41)). 
Only the last two equations should be changed, specifically 
Hose (LLet40a) and (11.4 ]). 

Eq. (II.40a) was a constraining equation which ensured 
the confinement of the undesired roots to the left of the 


Soe line for the case of two real undesired roots. [In 


Toure 
this problem it takes the form 


a- (a - by)’ 20. (IIT .42) 


If the undesired roots are complex?, the equation 
= Ee fey (Git .43) 


serves the same purpose. 

These two constraining equations ensure that the magni- 
tudes of the real parts of the undesired roots are greater 
than the magnitudes of the real parts of the desired roots. 


The performance index takes the simple form 


So: GC ae (II .44) 


because it is only necessary to maximize the real parts of 


Ehe desired@roots. 


1 For this problem all roots are undesirable, but just 


for purposes of terminology the names used in Example 2 
are retained. The roots of the polynomial P,g are still 
called "the undesired roots" and the roots of Pg are cal- 
led "the desired roots". 


Be, 


The results of the minimization are shown in Table 
II.3, and the time response to the unit step input is shown 


in Fig. I1.10. The desired response for Example 2 is also 


plotted in the same figure for easy comparison. 










UNDESIRED ROOTS DESIRED ROOTS SYSTEM PARAMETERS 





or | * * 
part jeralie! * Ko KF 
1000.} .0625 


TABLE II.3 


J. EXAMPLE 4 

In some design problems it may be necessary to insert 
a lag compensator which has a dipole (a pole and a Zero) 
very close to the origin. In this case some undesired 
root(s) may be located closer to the origin than the de- 
sired roots. Problems of this type may be solved by using 
the normal step-by-step procedure given in Example J..5 Spree 
this is fone an acceptable answer may be obtained by re- 
laxing the limiting value for the undesired roots (05 ,) by 
a considerable amount toward the origin (as in step 9 of 
Example 1). If the designer can gain enough insight to 
the problem (probably at step 2 of the normal step-by-step 
procedure) he may directly select Or, = O (to ensure the 
stability of the system) and save some computer time by 


avoiding the relaxation steps for Tre 
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| fe 


lee 
5 2 3 4. 2 


Fig. II.10. EXAMPLE 3. The fastest time response (A) and 
the desired response (B) for 
EXAMPLE 2, 


DUTPUT ——> 
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In another approach the designer may identify and con- 
fine some of the undesired roots to a close vicinity of jG 
Origin while keeping the other undesired roots in a remote 
region of the s-plane. In the design example given below 
this method is used just to explain a variation from the 
main method given in Example l. 

A third-order system (shown in Fig. II.11l(a)) is given. 
Steady-state performance requirements specify that the sys- 
tem gain (K) must be greater than 425 (the stability Lama € 
of the uncompensated system). The step-by-step solution is 


given below; paragraph numbers indicate the corresponding 


steps in Example l. 





(a) Uncompensated system (b) Compensated system 


Fig.II.1l. EXAMPLE 4 


1. Dynamic performance specifications require a pair 
of dominant roots located at O 5 = l, and w, =8e- 2 

2. - 3. The designer decides on a cascade lag com- 
pensator by sketching the root locus of the uncompensated 


system. The compensated system is shown in Fig. II.11(b). 


A2 


4. The transmission function of the compensated system 


is 


1 
s°+(12+p)s~+(36+12p) s°+(36p+K/k) s+Kp 


and the characteristic equation is 


P_» (@.8)= Ss" +(L24p) s°+(36+12p) s+ (36p+K/k) stKp —ag 


(II.46) 
where 
eo (et hes 
5. The quadratic for the desired root is 
P,(o,W,s) = ae + 20 ste ( dee4s7) 


where 


ch ot tw. 


For this problem the characteristic equation of the 
system cannot have more than two complex roots; since the 
- designer is interested in the part of the root locus where 
the desired roots are complex, then the undesired roots 
must be real. For this reason the polynomial representing 
the undesired roots can be written as 


Pg(Ays) = (sta)-(stb) = s*+ (atb)s + ab (11.48) 


The characteristic equation is factored by using id@- 


ELey (Lio 


P = Pp_p (EL s)) 
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¢*4+(12+p) s-+ (6412 5)sc (36p+K/k)s + Kp = (II.49) 
s*+(a+b+20) s°+(abt+e+2a0+2be) s-4+(2abo+actbe) s + abc, 


and 0, p, and b are chosen as the variables in the minimiza- 


tion process; in other words, 
T 
Z= (0 p b) (IT.50) 


The constraining equations can be written (from Eq. 


(II.49)) as 


a =o2)4+ p =—SeeeZo 
c = 36 + 12p - 20(atb) - ab 

( Tien) 
K = abc/p 


k = K/(-36p + 2abo + c(atb)). 


6. All of the variables must be greater than or 


equal to zero, and the other constraining inequalities are 


K 2 425 

[=k cS 26 (LESS2 ) 
b> Ore] limit for the remote undesired root 

Oe =" ag = 7a limits for the close undesired root. 


1t-2 
7. The performance index 


Z 


J = (1.-c) 7 + (2.-W) (Tig) 


is used for minimization by SUBROUTINE BOXPLX. 
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9. Some suitable logic to relax the boundaries was 
inserted in the main program. The remote undesired root 
Ctace——— os eonrmumnecad te the Jere of Ehe = = “F144 line 
and the close undesired root (at s = -a) is confined be- 
tween the s = O and s = “Ti 45 lines, where 1 E-1 = ie eval 
T1492 = 1.5. The first limiting line would be shifted to 
the right and the second one shifted to the left, but no 
relaxation is done because the minimum is not on the enclos- 
ing boundaries. The result of the minimization is tabulated 
as Table I1I.4. 

10. The compensated system's response (to a unit step 
input) and the desired response are shown in Fig. II1.12. 


PANO DeatncuGesd’ Hod roOte are Mocated exactly at tie ™Sclec] 


ted locations, they are obviously not dominant becasue the 


two time responses are not close to each other. 













LIMITING VALUES 
FOR THE UNDE- 
SIRED ROOTS 





UNDESTRED ROOTS | DESIRED ROOTS} SYSTEM 


PARAMETERS 







| -0, 


or fee 1 Close | Remote 
E." 2 rere 
TABLE I1I.4 
K, CONCLUDING REMARKS 


The method "Optimization for the Best Root Locations In 
the s-Domain" provides fast, direct and efficient solution 


to the linear system design problems of classical control. 
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EXAMPLE 4. Time response (A) of the compensated 


system with a lag compensator and 
the desired response (B) to a unit 
step input. 


46 


Systems with more than two free parameters can be handled 
easily. 

For high-order systems, algebraic manipulations neces- 
sary to factor the characteristic equation may get too in- 
volved. If the number of free parameters is very small 
compared to the order of the system, these algebraic mani- 


pulations may become prohibitively time consuming. 


A7 


IIIT. OPTIMIZATION FOR THE BEST RESPONSE IN CHE 
TIME-—DOMAIN -- APPLICATIONS TO LINEAR 


SYSTEM DESIGN PROBLEMS 


A. GENERAL 

Most of the methods for the design of dynamic systems 
are inherently indirect, s-domain methods, for example, 
try to adjust the root locations of the system character- 
istic equation to obtain a suitable transient response. 
For a high-order system it may be very difficult to decide 
on a set of suitable root locations which defines a desir- 
able transient response; for this reason, the final step 
of an s-domain design procedure is usually the simulation 
of the system, to check if the new root locations yield 
an J éceneanne transient ssesponse: 

The method presented in this chapter eliminates the 
intermediate steps and directly tries to shape the tran- 
sient response of the system. This method is called "Op- 
Pin teacon for the Best Response in the Time-Domain". 

The first step of the method is selecting a"MODEL 
RESPONSE"+ by interpreting the dynamic performance speci- 
Licatlons we laen ene: seace equations of the system are 


integrated with respect to time with a given input 


1 
The word "response" implies "time response" even when it 
is used alone. 
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function, and the output of the system is obtained and 
sampled at suitable time intervals. These samples are com- 
pared with the desired values of the response (defined by 
the model response). A performance index, which measures 
the deviation of the sample values from the desired values 
(from zero time to the expected end of the settling period 
of the system) is defined. Minimization of this performance 
index with respect to the free system parameters yields the 
best time response of the system, for the given specifica- 
tions and the system configuration, 

In this chapter, the simple philosophy of the method 
is described and the method is applied to the linear sys- 
tem design problems which were solved in Chapter II by 
using the method "Optimization for the Best Root Locations 


in the s-Domain",. 


B. SELECTING THE MODEL RESPONSE 

Selection of the model response is determined by the 
dynamic performance specifications. If the system is to be 
operated always with a specific input signal, then this 
time function (input) can also be used as’ the model re- 
sponse for the comparison. A step input (test signal) is 
very common in dynamic system design procedures. The same 
step function may also be used as the model response; if 
this is done, however, it may be necessary to insert some 
constraints on the states of the system to fulfill other 
dynamic performance specifications, such as peak overshoot, 


etc. (see Examples 7 and 9). 
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If the dynamic performance specifications are too 
tight (that is, if expected rise time, peak overshoot, set- 
tling time, frequency of oscillations are given within some 
small tolerance limits) a second-order model system may be 
used to represent the desired model response easily and 
efficiently. The output of the model system is obtained 
by integrating the state equations of the model. Using a 
second-order system's output as the model response is of 
course, equivalent to the dominant-roots concept in the 


s-domain. 


Cx CONSTRAINTS 

The constraints on the free system parameters are 
usually linear or can be approximated by linear inequali- 
ties -- this simplifies the computer programming part of 
the design procedure. 

There may be some constraints on the state of the sys- 
tem due to inherent nonlinear characteristics of the system 
components. These constraints can be represented by linear 
Or nonlinear inequalities and they do not introduce much 
difficulty to the design procedure (see Chapter IV). 

To obtain some design characteristics, additional con- 
straints on the states of the system may be introduced by 
the designer. For example, in one design procedure the 


peak overshoot of the system can be kept below a certain 


level by inserting suitable constraints on the output 
state. In this thesis, penalty functions are employed ze 


handle this type of constraint (see Examples 7 and 9). 
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De DEVELOPMENT OF THE METHOD 
The state equations of an n="_order, Single~input sys- 


tem can be written in the form 


x(t) = Ax(t) + br(t) Gait .1) 
Veehte mx (tlic athe state m-vector and the "dot" over a vari- 
able is used to indicate the time derivative (d/dt). Aes 
ann xX n matrix, r(t) is the scalar time function represent- 


mig the maput Co the system and b is an n-vector. 


The output equation is 
a 
c(t) =a x(t) + ne re) (rit. 2) 


where c(t) Po edio sOluS ibn Ot Joes) Son aS onl vec lok 
and ho is a scalar. 

Assume that the design is to be carried out by using 
a second-order model response. The state equations for the 


second-order model system can be written as 


O ih 0 
es) -y(t) + - r(t) (ieee 3 ) 
fd 2? P Fad 2 
~w" --20 w Ws 
na dn na 


— — 4 


lohis step may take different forms depending on the choice 
of the model response. If the input time function (r(t)) 
is to be used as the model response, for example, no extra 
effort is necessary to generate the model response. 


Sal 


where y(t) is the state vector of the model system, €. re- 


d 


presents the desired damping ratio and we represents the 
; 4 
desired natural frequency of the model system. The output 


of the model is 
c(t) = y, (t) P (Tli343 


A performance index which is suitable for many design 
problems is 


Ct f 
: 2 
s(@ = [ (eg (t). - e,(t)) ae , (IZI.5) 
s 

where @ 1s "an N-dimensional vector representing the free 
system parameters, and: c(t) represents the model response 
(for the case under discussion 1€ is the outpwe of@en 
second-order model). ot, (the final time of the integeas 


tion) can be taken as the settling time of the second- 


order model; it can be calculated by using the equation 
— we, a 
ia =e SS (III.6) 
d ns d 


where O54 is the real part of the desired root locations. 
If another model response is used, the final time can 
be guessed by using the performance specifications -- it 
should be long enough to cover at least one period of the 
lowest frequency of the transients which occur during the 


minimization process. The value of the performance index 


SZ 


depends on tee but once it is chosen (long enough) and 
fixed, it does not affect the optimum parameter values ob- 
tained at the end of the minimization. If ty 1s chosen 
longer thaniis necessary, extra computer time may be spent. 
For this reason, if the designer does not have any idea 
about the dynamic behavior of the system, he should simulate 
the system by using the starting values of the free system 
parameters as the first step of the design procedure. 
The performance index in the sampled form is written 

as 

oie 

) (ce (KP) = c (Kr))? (III.7) 
k=0 


J (&) 


where ne 1S a poSitive integer defined as 


_ final time of the calculations _ Xs 
i 


a Sampling period ° (ide) 


When J is minimized with respect to the free system 
parameters (&) in the bounded region defined by the con- 
straints (see Section III.C), the optimum parameter values 
(a*) are found. 

The form of the performance index can be modified to 
fit the needs of the specific design problem; variation 
from the basic form (Eqs. (III.5) and (III.7) are employed 
in the examples given in the later parts of the thesis 


(see particularly Examples 7 and 9). 
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Ky | EXAMPLE 5 

The linear system design problem solved as Example 1 
in Chapter II is reconsidered here. This time, a solution 
is obtained by using the method "Optimization for the Best 
Response in the Time-Domain". The system (shown in Fig. 
II.3(a)) and all specifications are kept the same to pro- 
vide a basis for comparison. A step-by-step solution is 
given below: 

1. The Model Response 

The desired dynamic response was represented by 

two complex roots, in Example 1; the desired dynamic re- 


sponse (the model response) is the output of a second-order 


model system defined by dominant roots located at a5 OS, 
w,= +t 0.5. The state equations of the second-order model 
are written from Eq. (III.3) 
x(t) = x, (t) 
2 (ITI.10) 


where x, (t) = y, (t) and x (t) = y(t) are used for program- 
ming convenience. A unit step input is assumed (r(t)= fl(t)). 


The model response is 


c(t) = x, (t) : (Titeaee) 


2. Type of Compensation 


The same cascade compensator (shown in Fig.II.4(a)) 


is used for compensation. 
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3. State Equations of the System 
The state equations of the system are written from 
the transmission function (Eq. (II.11) in the rational ca- 


nonical form 


x, (t) =X. (t) 


Xo (t) = x, (t) 


(ITE, ey 

x, (t) = x, (t) + h, (CE) 

x, (t) = “box, (t) - b,x, (t) - box, (t)- b,x, (t) + h, x(t) 
where 

bo = Kp b, =2+p 

b, =p + K/k h, = K/k (gett . 13) 

bo = | -=2p hy = bo ~ h,b, : 
The output of the system is 

c(t) - x, (t) A (III.14) 


4. The Performance Index 
The performance index defined in Eq. (III.5) is 
suitable for this design problem, hence 


a 


2 > _ 2 
J = x(t) = q (x, (t) x, (t)) Gey (TIT. 15a 


And ree Eime derivative is 


X(t) = (x, (t) - x_(t))* (Tir 1G 


The final time for the integration is taken as the 


settling time of the model response 


ty = 4/0 5 = 8 seconds. (Ilia 


5 BeCeonstraints 


The only constraining inequalities are 


O54 Syke 0. oe 


62a Manama Zataon 
The performance index (x  (t) ) is minimized with 

respect to the free system parameters (@)2 A gradicems 
method was used for the minimization (COMPUTER PROGRAM II, 
given at the end of the thesis is typical for this method). 
At each step of the minimization procedure, Eqs. (III.10), 
(III.12) and (III.16) are integrated from t = 0 to t. by 
using a fourth-order Runge-Kutta method~. The optimum 
parameter values are k* = 0.1, p* = 9.84, and K* = 0.4945. 
The optimum time response of the system and the desired 
response are shown in Fig. III.1l. 

The system's response shown in this figure (and its 
closeness to the desired response) should be compared with 


the response shown in Fig. I1.7. 


* SUBROUTINE RKLDEQ, N.P.G.S. Computer Facitlity Program 
Library. 
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Fig. III.1. EXAMPLE 5. Compensated system's response (A) 


and the desired response (B) to a 
unit step input. 
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Fe, EXAMPLE 6 

Example 2 of Chapter If is "re—-solved bys vcmaceene 
method "Optimization for the Best Response in the Time- 
Domain". The system is shown in Fig. II.8. 

In the solution given below, paragraph numbers indic- 
ate the corresponding steps in Example 5. 

1. The model response is obtained from the output of 
a second-order model system which is defined by the state 


equations 


X.(t) = x, (t) 


(III.19) 
xe (t) = 5(r(t) - x, (t) - 2x,(t) , 
and the model response is 
c ft) = x. (t) ; (ITi. 2) 
A unit step input is used in the design, thus 
r(t) = f(t) . (lee) 


2. The same tachometer and acceleration feedbacks are 
used to compensate the system. 

3. The state equations of the system are written 
directly from Fig. II.8(b). 


x, (t) = 35 (0) 


Xe UE) tf 
ae (TIT see 


X(t) = x, (t) 


x, (t) = [r(t)-x, (t)-K,x, (t) -K x, (t)] K-[50x, (t)+65x, (t)+16x,(t)1 
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The output of the system is 
c(t) = x, (t) : ( Lieg2i3) 


4. The time derivative of the performance index is 


_ I 
x7 (t) = (x, (t) - x, (t) ) : (IIL .24) 
The final time for the integration is 


t = 4/0 4 =a le eoee ones. (TEL .25) 


5. The constraints on the free variables are 


< 
Os KL 1O , 


(LIL. 269 
Os Ko = 10", 
and the system gain is set to its desired value (K = 1000). 
6. The performance index defined by the integral of 
Eq. (III.24) is minimized by using a gradient method; the 
optimum parameter values are Ks = 0.142, Ky = 0.4, and the 
optimum time response of the system (for a unit step input) 
is plotted in Fig. III.2. In the same figure the desired 


response is also shown. The result of s-domain synthesis 


procedure was given in Fig. II.9. 


Gas EXAMPLE 7 
The linear system design problem solved in Examples 
2, 3, and 6 is re-considered here. As in Example 3 of 


Chapter II dynamic performance specifications require "the 


De, 
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Fig. III.2. EXAMPLE 6. Compensated system's response (A) and 
the desired response (B) to a unit step 
amar u E < 
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fastest response” to a unit step input. This problem is 
solved three times by using slightly different techniques 
to demonstrate the effects SEepene dlrtevent forms Of per 
formance indices. 

So@fution 1 


1. .The unit step input is used as the model response: 
eres age) ee) ae (DDI. 2 


Other parts of the solution are exactly the same as 
the one given in the previous example except that the per- 


formance index now takes the form 


Sale) (L.=x, (t))? (III. 28) 


2. The new performance index te minimized with respect 
to the free system parameters (K, and KY) ant the optimum 
parameter values are K* = 020296, KE = 0.4546. The optimum 
response is plotted in Fig. III.3 . In the same figure the 
desired response of Example 6 is also shown, for comparison. 

Pouucion, 2 

The optimum parameter values in Solution 1 yield a 
faster response (shorter rise time and shorter settling 
time) than obtained in Example 6, but the response is some- 
What oscallatery and the settling time is about four sec- 
conds. If dynamic performance specifications require a less 
oscillatory output and a settling time less than, for exam- 


ple, 2.5 seconds, a small change in the design procedure 


yields the desired dynamic response. 
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(for EXAMPLE 6) 
(SOLUTION 1). 
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The fastest time response (A) and the 
to a unit step input. 


desired response (B) 


fli 
EXAMPLE 7. 


O 
5 ip eee 


Fuge 


For this application the performance index takes the 


f£omm 


EE 


J = x2 (t) = f w, (1.-x, (t)) “at ’ (TTT... 29) 


where t. = 4 seconds as before and Wy is the penalty func- 


tion defined as 


a moe to — 265 


“a = (LIL «3@) 
1000 GeEOG ee. 242.508 


This penalty function assures a settling time less than or 
equal to 2.5 seconds. Minimization of Eq. (III.29) yields 
the optimum parameter values; Ke = 0.047, Ke = 0.441. 
The output of the compensated system with optimum parameter 
values is shown in Fig. III.4. 

Solmeaon. 3 

A somewhat different approach to the same design prob- 
lem is considered. The differences from the previous 
solutions are given below. 

The system gain is used as the third free parameter 
(which is constrained in the region 10° SK < Oe Ne 

The performance index 


fe 
wee. 7 
_— x(t) = J | w, (L.-x, (t)) - x(t) | dt (TET 23) 


is used for the design. In this approach the time deriva- 


tive of the output is maximized to obtain the fastest 
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1p Be 1 





TII.4. EXAMPLE 7. 











The fastest time response (A) (with a 
constraint on the settling time) and 
the desired response (B) (for EXAM- 
PLE 6) to a unit step inpae. 
(SOLUTION 2) 
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response. The first term (w, (1. =x, (t))*) is used to satis- 
fy another (assumed) performance specification which requires 
the peak overshoot be less than or equal to 40%; for this 


purpose the penalty function 


1 for | ee (1) |< 14 


1 for all t in be, t 1 
100 fOr |somi(et2) ss i 4 


(rrr ap 
is used, 
Minimization of the performance index (Eq. (II1I.31)) 
yields the optimum parameter values K* = 1012.6, K* = 0.07, 


K, = 0.296. The optimum response is shown in Fig. III.5. 


re. EXAMPLE 8 

In the Example 4 of Chapter II a third-order linear 
system was compensated for the given dominant root loca- 
tions. The use of a lag compensator (in the s-domain) was 
demonstrated in Example 4. The same example is re-solved 
here to compare the results of the two methods. 

The system is shown in Fig. II.1ll(a). Paragraph 
numbers are used to indicate the corresponding steps in 
Example 5. 

1. The desired root locations in Example 4 (O5 = l, 
w= + 2) correspond to a second-order system which is 


defined by the state equations 


X(t) = x, (t) 
(TIT. 32) 
X(t) = Seer) — X(t) ) = 2x- (t) : 
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i] Ge ees 


EXAMPLE. 






a a eae ty Salles ae 


The fastest response (A) by maximiZing 
the derivative term, and the desired 
response (B) (for EXAMPLE 6) 

(SOLUTION 3). 


66 


The model response is 

cbt) = ge (lo . (TEL. sy 
The unit step input is used in the design, i.e. 

r(t) = f(t). (IIT. 34) 


2. - 3. A cascade compensator is used (as shown in 
Fig. II.11(b)) to obtain the desired response. The state 
equations of the system are written from the system's 
transmission function (Eq. (I1I.45)) in the rational canon- 


ical form 


x(t) = Ax(t) + br(t) (Gat . 3.5) 
where 
ee 
ii 0 0 
0 
A = 0 0 0 il 4 (TIIT.36) 
-a ~-a -a -—a 
Hf 2 3 
° 2 
Eee O00 her. amd 
~ SY 4 4 
a+ Kp a= Ve Se 72 
( Leni) 
a, = 36p+K/k h.= K/k 


fe) 
Il 


SG LZ p h,= aot h,a, 4 
The output of the system is 


c(t) = x, (t) : 


67 


4. The time derivative of the performance index is 

x(t) = (x. (t) - x.(t))? ? (ITI.38) 
i) 1 5 : 

5. The constraints on the free system variables are 


425 = K s 1000 
O.l1 =k <= 10 (lll eS 


O<p FelOOm(arbatrany) . 


6. The performance index defined by the integral of 
Eq. (III.38) is minimized with respect to he aes system 
parameters and the optimum parameter values are K* = 554, 
k* = 7.84, and p* = 0.0116. The optimum response to a 


unit step input and the model response are shown in Fig. 


PiieG, 


Le CONCLUDING REMARKS 

In this chapter, the method "Optimization for the Best 
Response in the Time-Domain" has been developed and applied 
to several linear system design problems. The same exam- 
ples, given in Chapter II for the s-domain design method, 
have been reconsidered here to make comparison of the two 
design methods possible. From this comparison and from 
the general development of the method, it may be concluded 
that this method has many advantages over the present dynam- 
ike. system design procedures. The most important advantages 


are listed below: 
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Fig. III.6. EXAMPLE 8. 


TIME——~ 


Time response (A) of the compensated 
system with a lag compensator and the 


desired response (B), 
Io uit. 


69 


to a unit step 


1. This method is applicable to high-order, multi- 
parameter systems. Any dynamic system which can be repre- 
sented by state equations can be designed by this method; 
the order of the system and the number of free system para- 
meters can be very high without introducing any difEfieumiie 
in the design procedure. 

2. The work required to design a high-order system is 
less than with other methods; in fact, if a computer pro- 
gram for the minimization process is available, the only 
required work is in obtaining the state equations of the 
system. 

3. This method directly shapes the time response of 
the system by adjusting the free parameters. If the op- 
timum response obtained by this method is not close enough 
to the desired (model) response, then the system configura- 
tion must be changed. The only indirect part of the method 
is that the designer must decide on the available system 
components and compensation techniques. 

The method "Optimization for the Best Response in the 
Time-Domain" enjoys the following advantages over the s- 
domain method developed in Chapter II. 

a. The basic idea behind this method is much sim- 
pler. 

b. Optimization for the best response in the time 
domain yields a better answer to the design problem because 
it is more direct and it does not employ any questionable 


assumptions such as dominancy of the desired roots. 
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c. This method requires less hand calculation and 
less programming effort. 

d. This method is not limited to the design of 
linear systems. Besides its simplicity this is the most 
important advantage of the method, because it provides a 
powerful computerized approach to the design of nonlinear 
and sampled-data systems and also, for example, systems with 


time delay. These applications are presented in Chapter IV. 


Gt 


IV. OPTIMIZATION FOR THE BEST RESPONSE IN THE 
TIME DOMAIN -- APPLICATIONS TO NONLINEAR 
AND SAMPLED-DATA SYSTEMS AND SYSTEMS 


WITH TIME DELAY 


A. GENERAL 

The method "Optimization for the Best Response in the 
Time-Domain" is applicable to a large variety of dynamic 
system synthesis problems. Any dynamic system which can 
be represented by state equations can be designed to yield 
a desirable dynamic and steady-state performance. 

In this chapter, applications of the method to non- 
linear and sampled data systems and to systems with time 
delay are considered. Three simple, but representative, 
examples are given. Application to these classes of sys- 
tems increases the value of the method, because direct and 
efficient design methods are rare in these areas of control 
theory. 

If the system is not linear, it may require some addi- 
tional logic (in the computer programming) to obtain the 
state values in the time interval of interest; but the 
main procedure (as presented in Chapter III) is directly 


applicable. 


Bis EXAMPLE 9 
A chemical process with an inherent time delay (T = 0.5 
seconas) is to be regulated. The system is shown in Fig, 


IV.1. The load changes are sudden and can be assumed to be 
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step functions. The only available free parameter is the 
amplifier gain (K.)- The desired dynamic and steady-state 
specifications will be defined later. 

For the given system configuration the system's re- 
sponse can be divided into two parts: 

1. Steady-State Accuracy. Since the given system is 
type-0O, there will always be a constant steady-state error 
for step disturbances. The value of the steady-state 


error (E_.), which depends on the amplifier gain (Ki), is 
EK) = 1/(1+4.8 K) Civ") 


2. Dynamic Accuracy. This is also a function of the 
amplifier gain (K.)- The dependence of the shape of the 
dynamic response (on K.) will be apparent at the end of the 
selution. 

In general, high gain values cause more oscillations, 
and less steady-state error (see Eq. (IV.1)); low gain 
values result ina less oscillatory response but more 
steady-state error. With only one free parameter, it is 
impossible to meet all of the specifications; therefore, 
the design method should seek the best compromise between 
the desired specifications. 

To make the initial conditions of the states all zero, 
it can be assumed that the reference signal is zero (u(t)#=0)), 
and that the system is in a steady-state condition. The 
system's behavior in case of sudden load changes will be in- 


vestigated. 


Ts 


REF. LEVEL LOAD 


U(s) 





MEASURING DEVICE 3 


Fig liv. 1. «BXAMPEEM 


thesstate equations of thewspysitemyare 


x, (t) aioe ttt —0). Se 
( Tee) 


-(9.6K_ x, (E-0.5) +2x, (t) +9x, (t)+r(t-0.5))/9 


x, (t) 


There are three terms affected by the 0.5 - second 
time delay. x, (t-0.5) represents the value of x, (t), O25 
seconds earlier during the process. The load disturbances 


are assumed to be unit step functions; therefore, 


fe for t < 0.5 


mi E=0e >) = (iVeen 
| t fOr” Cen go 
can be used. 
The output of the system is 
(IV.4) 


c(t) — x, (t) . 


It the performance specifications only require that 


the output of the system should be kept as close to the 
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reference as possible, and the transients should last not 
more than 30 seconds, then a performance index, which meas- 


ures the deviation of the output from the reference level 


Giron t = GP te t p= 30 seconds) can be defined as 
IE 
a 5 
o(K,) 2 x,(t) 2 | Came)" — uleyiiede®. Cis 
O 


Since the reference level was assumed to be zero, this per- 


formance index becomes 


t 


5 UG le [@ 


1° (t) ole (IV.6) 


The time derivative of the performance index is 


(8) = x,°(t) (IV.7) 


Minimization of this performance index with respect to 
K. yields the "best" response of the system within the 
limitations of the system configuration. 

Eqs. (IV.2) and (IV.7) are integrated with respect to 
time to calculate the value of the performance index at 
each step of minimization process. The integration step 
size is chosen as an integer submultiple of the time in the 
process (At = 0.05 is chosen). For the first ten integra- 


tion steps the values of r(k At - 0.5) and x, (k At - 0.5) 


eS 


are eee the value of the variable x, (k At) are stored 

(k is the integrating step counter). After the tenth in- 

tegration step (k > 10), r(k At - 0.5) is equated to unity 
and the values of x, ((k-10) At) are found and used from the 
Stored valuccmou x, (k ISSO 

The only difference from a standard linear system de- 
Sign problem is in the logic of expressing the existing 
time delay in the system, and as explained above, this does 
not introduce any complication to the solution. 

K* = 1.66 is found by the minimization procedure. The 
time response of the system to a unit step load change is 
shown in Fig. IV.2. Peak overshoot is about 33% and the 
steady-state error is 11.1%. 

If the desired response is to be less oscillatory, one 
can include the time derivative of the output in the perform- 


ance index; that is, 
<= x aie) + xX ate) (IV.8) 
S li ila A " 


or by substituting the value of x, (t) (£rom Eq. (IV.2)) one 


obtains 


2S a ts) : (xs) + r(t - 6@5)/3)° . (ive 


: Since the numerical integration is a discrete process, 
the arguments of the delayed variables, r(t-0.5) and 
x, (t-0.5) are represented in discrete form. 
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Fig. IV.2. EXAMPLE 9. Time response of the system to 
a unit step disturbance (J = xy). 


a 


eee eee ee ae ——— 


By minimizing this performance index K* — lee 1S Ob= 
tained (The time response is shown in Fig. IV.3). This 
time response is less oscillatory but suffers more steady- 
state error (12.2%). 

To meet the given specifications, the terms represent- 
ing the output and its time derivative (in the performance 
index) may be multiplied by suitable weighting factors. 

For example, K* = 1.12 and the time response shown in Fig. 


IV.4 are obtained by minimizing the performance index 
x(t) =w, x,7(t) + wy x,7(t) (IV.10) 
3 Ihe a : = 


if @i2 Wy = land Wo = 10. 


If the desired response should have the smallest cuq 
mulative rate of change (slowest varying transients in the 


time interval yuipts. J!) the performance index takes the form 


aE) = x(t) = Gale) = ett — Se) (IV.11) 


This performance index is the same as (IV.10), of course, 
when Ww, = O and w, = 1 is used. With this performance in- 
dex K* = 0.473 is obtained (for the time response see Fig. 
TVS oes 

Finally some applications may require that the system's 
output should: not deviate from the reference level more than 


a certain percentage for a given load change. To satisfy 


this requirement a constraint on the output is inserted. 


WAS 





— 


Fig. IV.3. EXAMPLE 9. Time response of the system to a unit 
step disturbance (J = x1° + x1°). 
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Fig. IV.4. EXAMPLE 9. Time response of the system to a 
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TIME 





30. 40. 


Fig. IV.5. EXAMPLE 9. Time response of the system to a 
unit step disturbance (J=x5). 
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For example, to ensure that the maximum deviation from the 
reference level is less than 31% of the load disturbance, 
the performance index defined in Eq. (IV.8) is used and the 


additional constraint 
ez Ce) | < 0.31 forall t in Hon Ee! (IV.12) 


1S imposed. Minimization” yields Ke = 1.85; the time re- 
sponse of the system (for a unit step load change) is shown 
mn bac 1 a6 | 

This series of solutions (see Table IV.1) shows tao 
a system can be optimized for any desired response (that is, 
the most suitable response within the limitations of the 
given system configuration can be obtained) by using the 
method "Optimization for the Best Response in the Time- 
Domain" -- the existence of one or more time delays does 
not introduce any difficulty or complication in the solu- 


tion process. 


The penalty function technique used in Example 7 was 
applied. 
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Fig. IV.6. EXAMPLE 9. Time response of the system to 


a unit step disturbance 
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TABLE IV.1. EXAMPLE 9 (SOLUTIONS) 


: EXAMPLE 10 


Q) 


A simple second-order sampled-data system (shown in 
aCe. 7 Gee) eens considered. The two adjustable design 
parameters (free system parameters) are the system gain (RK) 


and the sampling Qecod CT) 


1 The method "Optimization for the Best Response in the 
Time-Domain" yields an important advantage when it is ap- 
plied to the sampled-data systems -- the sampling period 
can be taken as one of the free variables without increas- 
ing the complexity of the solution. 


! 
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Fig. IV.7(a). EXAMPLE 10. The system. 
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Solution (without a Hold Cisse uae) 

1. The dynamic performance specifications* require 
that the system's response should be as close to the re- 
sponse of the second-order continuous system as possible. 
The state equations of the second-order continuous system 


“i 


are 


x(t) = x, (t) 
; ‘i (Ty Bey 
x, (t) “= 0.526(r(t) - x,(t)) - 0.630x,(t) 
and the model response is 
cb) = x (t) : (IV.14) 


The unit step input is used in the design, i.e. 
vel((e)) ge oad 2) (TVs) 


3. The state equations of the system are written from 


Fig. LVe¢ay 


x, (t) = x, (t) 
(IV.16) 
Xo (t) = K e, (t) Srey) ae 


Since the hold circuit is not used in this solweaen 


(see Fig. IV.7(a)). 


Paragraph numbers indicate the corresponding steps in 
Example 5. 


S6 


e, (t) - eae): (Tali?) 


and also 


0 if the sampling switch is open 


e, (t) 
r(t)-x, (t) during the sampling instant 
(IV.18) 

Logic must be included, in the computer program, to 
find the sampling instants, during the integration; in 
fact this is the only required addition to the computer 
program which is written for a standard linear system de- 
sign. 

It is assumed that the sampling switch is closed for 
the first time at t = 0, and opened after an integration 
step. In other words at the first step of the discrete 
numerical integration process e, (t) = r({t) - x, (t) is used, 
then r, (t) is set to zero. Numerical integration is con- 
tinued for "I - 1" steps with the e, (t) = 0 vakbue. I is an 


integer found from the equation 


_ Sampling Period = 
v Integration step size T/At (acs IED), 


A step counter (i) counts the integration steps starting 
from the sampling instant. i= I - 1 indicates another sam- 
pling instant: when this occurs 

ie = BOG 


(IV.20) 
e, (t) =e(t) - x, (t) 


BF 


are set by the computer program, and a new count is started, 
The sampling time 1s assumed to be one integration step size 
in this solution. If the sampling time (the period of time 
that the sampling switch stays closed) should be long com- 
pared to the integration step size then another counter may 
be added to indicate the opening instants of the sampling 
switch. The integration step size should be chosen as an 
integer submultiple of both the sampling period and sampling 
time. 


The output of the system is 


c. (t) = x, (t) . (IV e2a5) 


4. The time derivative of the performance index is 
a(t) = (x, (t) - x4 (t))7 (IV. 22) 
5 ih 3 ° . 


The final time of the integration is taken as the set- 


tling time of the continuous model, 


t = 4/0 5 = 12.7 seconds . (IV.23) 


5. Constraints on the free variables are 


O.L Se = ae 

(IV.24) 
eK s fiw 
6. The performance index defined by the integral of 


Eq. (IV.22) is minimized with respect to the free system 


parameters (T and K) and the optimum parameter values are 
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T™ = 0.6 ; kK" = Feo. The output of the system to a 


unit step input and the model response are shown in Fig. 
TV. 7a) 

Solution (With a Hold Circuit) 

If a zero-order hold circuit is inserted after the 
sampling switch the only necessary change is in the logic 
given in conjunction with the Eq. (IV.18). 

After a sampling instant, e, (t) is not equated to 
zero, but it is retained at the value of the last sampling 
for the remaining I - 1 steps of integration. [In this 
solution impulse sampling is used, i.e., the switch is not 
kept closed for 0.05 seconds but it is opened as soon as 
the sample is measured. 

Minimization of the same performance index (Eq. (IV.22)) 
yields T* = 1.05 and K* = 0.5875. The output of the sam- 
pled data system with zero-order hold and the model response 


are shown in Fig. IV.8. 


Dy EXAMPLE 11 

A designer Simulates the fourth-order linear system 
shown in Fig. IV.9 and obtains the time response to a unit 
step input shown in Fig. IV.10. This response is for a 
system gain of K = 1600. It is assumed that this dynamic 
response perfectly fits the particular application,except 
that the velocity coefficient must be doubled for steady- 
state accuracy. When the gain is doubled (K = 3200) however, 
the system moves to its stability limit, and the time re- 


sponse is quite oscillatory as shown in Fig. IV.10. 
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Fig. 


IV.8. EXAMPLE 





10. 


318 


Time response of the 
system (A) with a zero 
order hold circuit, 
and the desired 
response (B). 





R(s) 1600 (s + 8) C,(s) 


7" s(s+5) (s+10) (s+20) 7 


Fig. IV.9. EXAMPLE 11. The linear system. 






The real system has two inherent nonlinearities -- a 
saturating amplifier in the forward path and dead zone in 
the feedback path. The real system with the nonlinearities 
1s shown in Fig. TA aioe The designer decides to use two 
stages of linear cascade compensators to obtain a time re- 
Sponse as close to the response of the linear system as 
possible; in other words, the linear system shown in Fig. 
IV.9 (with K = 1600) will be used to generate "the model 
response" for the optimization process, 

The compensated system has four free parameters (see 


Fig. IV.12) -- k k, (the compensation ratios of the first 


- 2 
and second compensators) and Py Po (the pole locations of 
the compensators). 


The state equations of the ideal (linear) system are 
x(t) = x_(t) 
X_(t) = x(t) 


p (LVe25) 
X(t) = Xy oft) + K. (r(t) - x. (t)) 


I 


Kft) = -(1000x, (t) +350x, (t)+35x,, (t)+27K, (x(t) ~x, (t)) 


eal 


(A) 
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Fig.IV.10. EXAMPLE 11. Time response of the nonlinear system 
(A) (for K = 3200), and the time ree 
sponse of the linear system (B) (for 
K = 1600), (unit step input is used). 
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where the PLANT is defined by the transfert function 


c.(s) s+8 
Z (s) = Ss) = s(s+5) (st10) (s#20) 


(a) Real (Nonlinear) system 


Z(t) 





K = 3200, K = 1600, K =-1400 
x n 
(b) Ny - Block (c) N, - Block 


Fig. IV.11. EXAMPLE 11. The real system. 
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where K. = 1600 and r(t) = f(t) -- the unit step function. 


The model response is the output of the linear system 
c(t) = xo (t) : (CIN 2G) 


The state equations of the real (nonlinear) system 
with two cascade compensators (refer to Fig. IV.12) are 
written in the following steps: 


1. The plant is specified by 


x, (t) = x, (t) 


X5(t) = x,(t) 
Civic 27) 
X(t) = x, (t) + z(t) 
x, (t) a ~(1000x. (t) +350x,(t)+35x, (t) +27 Z,(t)), 
and the output of the system is 
c(t) = x, (t) . CIES) 


2. The Feedback Loop (N. - Block) 
The output of the dead zone block (see Fig. IV.11(€)) 


can be represented as 
O Wie PSC i 
Z(t) = x, (t) -K_ ai x, (t) > K. (IV.29) 


~ 
x, (t) +K, aE x, (t) K : 


So Wnentey 


The error is given by 


e(t) = r(t) - Z,(t) : CIV. 0) 


95 


4. Compensators 
The cascade compensators have transfer functions 


of the form 
G (s) = oo, ee 1 =¢,2 (IV seas) 
Cc, ‘og e a e 
at a 
in the time domain; they can be represented by the confi- 


guration shown in Fig. IV.12. For the first compensator 


fe (t) = -p, xXg(t) + py(1-1/k,) e(t) , (IV. 32) 


and its output is 


y,(t) = e(t)/k, + x, (t) . (IV. 33) 


Similarly, for the second compensator 
xa(t) = — pate. (Een enemy (IV. 34) 


and its output is 


y(t) = y,(t)/k, + x, (t) . (IV.35) 


Se Ny = Brock 


The output of the saturating amplifier is a non- 
linear function of the magnitude of its input (see Fig. 


IV.11(b)); its output can be represented as 


> 
K, af Z, (t) Ky 


Z, (t) = Ko ae Z(t) es (IV.36) 


_K y,(t) otherwise 


YS 


A performance index, which measures the difference of 
the outputs of the nonlinear and ideal systems, for a given 


magnitude of a step input (R), is defined as 


G 


fn 


T(R,@) = x,,(t) = f ogee - x, (t)) “at : (IV. 37) 
and minimized with respect to the free system parameters 
(Y% = (k, ko Py p5)”) and the optimum parameter values (a) 
are Louie), 

Since the real system contains nonlinearities which 
are amplitude dependent, the performance index and the op- 
timum parameter values are functions of the magnitude of 
the step input (R). For various amplitudes of step inputs 
the optimum parameter values are tabulated in Table IV.2. 
The time responses of the ideal and real systems for the 
conditions given in the last three lines of Table IV.2 are 
shown in Figs. IV.13, IV.14 and IV.15. 

If it is known that the actual system is going to be 
operated usually with step inputs of a certain amplitude 
(say, R= 2.5), then the design can be completed by using 
the last line of Table IV.2; in this case the dynamic re- 
sponse of the system may not be satisfactory for step inputs 
much larger or smaller than the selected value (R = 2.5). 
Mo study Oe effect of the input magnitude, the pole loca- 
tions are fixed at the optimum value found for the step 


input with 2.5 magnitude, 


De=wlOe2764 = Constant , (IV. 38) 
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TABLE IV.2 


STEP INPUT 
MAGNITUDE 


OLS 
i. SO00s 
i 7 ioe 
1.7480 
We? 25 
1.6904 
1.2132 
dee2 LOO 
d Fee 0) ial 
l-1ses 
et s02 
dL cules 
aati ss a 
Leo 
dls) 7): 
treo LOZ 
1.3543 
L7-Se, 


a al 
‘Say 
Us 
0.4 
Or 
0.6 
Ory 
amare: 
aro 
les.0 
le 
Ze) 
ae, 
Bi) 
oeres 
4.0 
4.5 
Dsl) 





TABLE IV.3 


98 





Fig. IV.13. EXAMPLE 11. Compensated (A) and desired (B) sys- 
tems' responses (r(t) = f(t)). 


NS, 








Fig. IV.14. EXAMPLE 11. Compensated (A) and the desired (B) 
systems' responses (r(t) = 2 x fl(t)). 
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Fig. IV.15. EXAMPLE 11. Compensated (A) and the desired (B) 
systems' responses (r(t) = 2.5 x @(t)). 


and the optimum compensation ratios (Ek) are found and 
tabulated (see Table IV.3) for various step input magnitu- 
des (R). 

A graph showing the relationship between the step in- 
put magnitudes and the compensation ratios (R vs. (ky and 
k5)) is plotted (see Fig. IV.16). This curve can be divid- 
ed into three regions: 

(l) Wel Seka 026 
(2) Oke AS Rasa (IV.39) 


(3) 9325 = Ree 


where R is the magnitude of the step input. The flat sec- 
ond region indicates that if a suitable value for the com- 
pensation ratios is chosen, the system has satisfactory 
output responses for the input magnitudes in the second 


region. For example, 


k, = ky 


Pi P2 


ie elae 
(IV.40) 


lI 


10.2764 


may be chosen as a (fixed) set of design parameters for the 
input magnitudes in Region 2. If this range covers most of 
the operational input values of the system and some degrada- 


tion from the desired response can be tolerated* for tie 


a 
Since the system has a stable limit cycle, growing oscil- 
lations do not occur for any input magnitude. 


loa 


Big. 
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IV.16. EXAMPLE ll. 
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Input magnitude (R) v.s. Optimum 
compensation ratios (K* and KS 
for pole locations are’ fixed 


at Py = Py = 10.264. 
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other input magnitudes, this design is acceptable; other- 
wise another compensation scheme (probably adaptive or non- 


linean) must be Considered, 


Ee CONCLUDING REMARKS 

The method "Optimization for the Best Response in the 
Time-Domain" is directly applicable to the design of sampled 
data systems and systems with one or more nonlinear elements 


or time delays. 
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V. A TIME-DOMAIN STABILITY ANALYSIS METHOD 


USING NUMERICAL TECHNIQUES 


Ae GENERAL 

A numerical method for the stability analysis of feed- 
back control systems is introduced in this chapter. This 
numerical method is applicable to a large variety of feed- 
back control systems, but it may be valuable especially 
for the nonlinear and sampled-data systems and systems with 
time delay, since general and practical methods are rare in 
these fields. 

If the stability analysis is carried out for a system 
which has been designed by using the method "Optimization 
for the Best Response in the Time-Domain", then only minor 
changes in the computer program yield the stability limit 
of the system; for this reason, this method is especially 


convenient for these types of design problems. 


Be GENERAL PHILOSOPHY 

The theory behind the numerical method is quite simple. 
It again depends on the idea of "shaping the transient re- 
sponse". The state equations of the system are integrated 
with respect to time by using a numerical method, the out- 
put of the system is observed during the integration (from 


zero time to final time - ty), and the maximum and minimum 


LOS 


points of the output are found. A performance index, to 
measure the sum of the differences of all the maxima and 


the differences of all the minima is defined 


a ee j-1,2 on j-1,2 
max “max es (Xin xmin ) (V.1) 


J(K) = 
where K is assumed to be the only free variable, the super- 
script j is the counter for the maxima and minima, and n 
is the number of maxima and minima found in the time inter- 
val of interest. 

Minimization of this performance index with respect to 
the free variable (K) yields the value of K*, which puts 
the system at its stability limit, because with the Ke value, 
the output of the system becomes constant-magnitude oscil- 
lations. 

For systems with two variables, a stability curve can 
be obtained by holding one of the variables at a constant 
value and applying the above procedure; at the end of the 
minimization, one point is found on the stability curve (in 
the system's parameter plane). By assigning other values 
to the fixed variable and minimizing the performance index 
with respect to the free variable, many points of the sta- 
bility curve can be obtained and plotted. 

The same idea can be extended to systems with more 
than two variables and stability surfaces or stability 


hypersurfaces can be obtained, but this may require too 
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much computer time, and the practical use of stability 
hypersurfaces in n-dimensional Euclidean spaces is ques- 
tionable. For practical design and analysis problems, the 
intersection of the stability hypersurface with a two-di- 
mensional plane may yield the desired insight, and it can 
easily be obtained. 

The terms in the performance index (values of the suc- 
cessive maxima and minima) are found by adding simple logic 
to the subroutine which calculates the value of the per- 
formance index. The maximum and minimum points of the re- 
sponse can be observed by different methods -- the simplest 
logic is discussed here. At the end of each integration 
step (of the system's state equations) the slope of the 
system's output is checked. If its sign is altered during 
the integration step, or if it takes the value zero, a max- 
imum or a minimum has occurred. For a continuous output 
response every maximum must be followed by a minimum and 
Vice versa; for this reason the performance index defined 


in Eq. (IV.1) may be altered to the simpler form 


n 
g(K) = ) xd = x72)? , (V.2) 
j=3 


where x represents both the maxima and the minima, and n 
is the sum of the number of maximum and minimum points 


observed in the time interval of interest. 


OF 


To demonstrate the application of the numerical sta- 


bility analysis method three examples are given. 


oF EXAMPLE 12 

The stability limit of the chemical process with in- 
herent time delay is to be found. The system (Fig.IV.1) 
and its state equations (IV.2) are given in Example 9 and 
will not be repeated here. 

A performance index as defined in Eq. (V.2) is mini- 
mized with respect to the only free system parameter (Ko - 
the amplifier gain) and K* = 2.125 is found for the sta- 
bility limit. To check this result, the system is simulat- 
ed with this gain value and the response shown in Fig. V.1 


is obtained. 


DD EXAMPLE 13 

A sampled-data system (with a zero-order hold) ‘was 
optimally compensated for a desired time response by using 
two free system parameters -- T (the sampling period) and 
K (the system gain) (see Example 10). The system is shown 
in Pee y 27 ae 

The sampling period can be varied from 0.1 second to 
1.5 seconds and the gain can be varied from zero to 15. 

The value of the sampling period is varied in steps; 
it is first set to 0.1 then increased by 0.1 steps (up to 
1.5). At each value of the sampling period, a performance 


index similar to the one defined in Eq. (V.2) is minimized 
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2 ek 2 5) Bs 


a 


Digital simulation of the system 
at its stability limit (K 


Fig. V.1. EXAMPLE 12. 


LO}, 


with respect to K and the stability curve in the T = K 
plane is plotted (see Fig. V.2). This curve simply gives 
the maximum allowable gain for a selected value of the 
sampling period, or vice versa. Several points on the sta- 
bility curve are checked by digital simulation; for example, 
T=0.2 , K = 10.34 are read from Fig. V.2 and with these 
parameter values the system is simulated, and the time re- 
sponse shown in Fig. V.3 is obtained. Similarly, for 

T= 1.5 and K = 1.74, the time response shown in Fig. V.4 


1s obtained, 


Big EXAMPLE 14 

The stability limit of the system, shown in Fig.V.5(a) 
is to be studied. The uncompensated linear system is un- 
stable when K > 193.23 (see’ Example 2). The system is com- 
pensated by using a tachometer, and Ky = 0.4 is selected 
for the desired response. The linear compensated system 
(with K, = 0.4) is unstable when K > 1020. 

The tachometer used in the system has a nonlinear 
characteristic and it saturates when its output reaches the 
saturation level (2). In other words, the output of the 


tachometer is a function of the magnitude of its input and 


can be represented as (see Fig. V.5(b)) 


8 e > 
a este Kc. (t) aa 
ig Cl) -Z_ abag Kc. (t) < -Z) (V.3) 
K,¢_ (t) otherwise 
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Fig. V.2. EXAMPLE 13. The stability curve in the "T - K" 
plane. 
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Time response of the system at 
one of its stability limits 
= 0.2, K = 10.34) 
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Fig. V.4. EXAMPLE 13. Time response of the system at 
one of its stability limits 
M = 1.5 , K = 1.74). 
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(a) System (the plant is same as shown in Fig. II.8(a)) 
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Fig.V.5. EXAMPLE 14. 
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For a given saturation level (Z,), the stability of 
the system depends on the value of the gain (K) and also on 
the magnitude of the test (input) signal (R). 

The state equations of the system, written directly 


from Fig. V.5(a), are 


seta acco (E) 


Xo (t) x (t) 


. (V.4) 
x,(t) = x, (t) 


x, (t) = K(r(t)- z(t)- x, (t)) - (50x, (t) +65x,(t)+16x, (t)), 
where the input is a step function with magnitude R, i.e. 
r({e = Rx aie. ey. 5) 
The output of the system is 
e_(t) =etpfe) (V.6) 
and its time derivative is 
eACo = xe) 2 (eh) 


For a given step input magnitude (R), the state equa- 
tions of the system are integrated by using a numerical 
method; after each integration step, the logic given in 
(V.3) is used to calculate z(t) (Z_ = 0.5 is used). The 
sign of x, (t) is also checked against its previous sign and 
the maximum and minimum points of the output (x, (t)) are 


found. The performance index (Eq. (V.2)) is minimized with 


ies 


respect to (K). The result of this minimization is the 
maximum allowable gain for a given input magnitude (R) and 
for the given saturation level (2 = 0.5). After a minimi- 
Zation is completed, another magnitude for the step input is 
chosen and the minimization is repeated. The maximum gain 
values (K) for each magnitude of the step input (R) are 
stored and the stability curve in the R - K plane (see Fig. 
V.6(a)) is obtained (for Zn = 0.5). When the input magnitu- 
de is small (R * 0.6), the operation of the system is linear 
and the stability limit is at K = 1020. As the input magni- 
tudes increase, the tachometer feedback becomes less effec- 
tive, due to saturation. For example, when the step input 
magnitude (R) is 100, the system is unstable for K > 196.07; 
for this reason the stability curve (shown in Fig. V.6(a)) 
starts from K = 1020 for small input magnitudes, and asymp- 
totically approaches the uncompensated system's stability 
limit (K = 193.23) as the input magnitudes increase. The 
same stability curve of Fig. V.6(a) is replotted to show 
the effect of magnitudes up to R = 100 (see Fig. V.6(b)). 

As a different problem, the input magnitude is fixed 
aie 

R = 2.5 = Constant , (V.8) 
and the effect of the saturation level on the stability 
limit 1sm@studied. 

The saturation level of the tachometer feedback (2) 


is kept constant, and the minimization is carried out with 
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Fig. V.6(a). EXAMPLE 14. The stability curve in the "R —- kK" 
plane. 
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_ Fig. V.6(b). EXAMPLE 14, The stability curve in the extended 
"R = K" plane. 


ibd: 


respect to K; when the value of K* for a given Ze is £ound, 
Z_, 18 changed to another value and the minimization is re- 
peated. The maximum K values for each saturation level on 


the tachometer feedback, in the range 
On. S72 Sore, | (V. 99 
m 


are found and the stability curve in the ZK plane (for 


ye) 
II 
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o. CONCLUDING REMARKS 

Two words of caution about the computer programming 
are appropriate here. The performance index defined in 
Eq. (V.2) has at least two minima; one of these occurs 
when the system is overdamped and the time response has no 
Critical points. In this situation both the value of the 
performance index and its gradient become zero, and this 
erroneously indicates a point on the stability curve. Some 
provision should be added to the computer program to avoid 
this erroneous result. Counting the number of maximum 
points in the time response, and finding new starting 
points if this number decreases below a certain predeter- 
mined constant, is one way of avoiding this difficulty. 

Secondly, the performance index defined in Eq. (V.2) 
measures the differences in successive peaks and successive 
low points of the time response; but it is also a function 
of the number of maxima and minima which occur in the time 


imterval of interest. This time interval may contain 
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Fig. V.7. EXAMPLE 14, The stability curve in the "Z - kK" 
plane. i 


120 


different numbers of maxima and minima for different para- 
meter values. To make the performance index independent of 
the number of maxima and minima, only a certain number of 
them may be included in the performance index, for example 
the third, fourth and fifth peaks and the low points of the 
time response may be measured and included in the perform- 
ance index. Averaging the performance index over the total 
number of maxima and minima may also overcome this difficul- 
EY « 

The two programming difficulties discussed above are 
not really serious and with a little care they can be over- 
come. The ibe al method used to investigate the stabil- 
ity of the feedback control systems can be very attractive 
for the stability analysis of nonlinear and sampled-data 
systems and systems with time delay. 

This method is an extension of the design method 
"Optimization for the Best Response in the Time-Domain". 

If a system is designed by uSing this method, the stabil- 
ity analysis can be carried out with little additional ef- 
fort, since only minor changes ae in the perform- 
ance index) convert the subroutine which was used in op- 
timal design to the one to be used in the stability analy- 
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VI. CONCLUSIONS 


Two synthesis methods for the design of multi-parameter 
dynamic systems have been presented. The powerful mathemat- 
ical tools of Optimal Control Theory have been applied to a 
variety of system design problems. 

The first method, "Optimization for the Best Root Loca- 
tions in the s-Domain" has been applied to linear system 
design problems. The second method, "Optimization for the 
Best Response in the Time-Domain" has been applied to the 
design of linear, nonlinear and sampled data systems and to 
a system with time delay. Eleven design examples have been 
considered; successful solutions of these examples may al- 
low one to conclude, "These design methods provide fast, 
direct and efficient solutions to practical system design 
problems, and they can compete with any existing design 
methods." 

The second method has been extended to a numerical sta- 
bility analysis method, and its effectiveness has been de- 
monstrated by three examples. 

The methods presented in this thesis provide computer- 
ized design procedures for a wide class of dynamic systems, 
and hopefully, they provide one of the necessary bridges to 
narrow the existing gap between the theory and practice of 


system control. 
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